Polaron to molecule transition in a strongly imbalanced Fermi gas 
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A single down spin Fermion with an attractive, zero range interaction with a Fermi sea of up-spin Fermions 
forms a polaronic quasiparticle. The associated quasiparticle weight vanishes beyond a critical strength of the 
attractive interaction, where a many-body bound state is formed. From a variational wavefunction in the molec- 
ular limit, we determine the critical value for the polaron to molecule transition. The value agrees well with 
the diagrammatic Monte Carlo results of Prokof 'ev and Svistunov and is consistent with recent rf-spectroscopy 
measurements of the quasiparticle weight by Schirotzek et. al. In addition, we calculate the contact coefficient 
of the strongly imbalanced gas, using the adiabatic theorem of Tan and discuss the implications of the polaron 
to molecule transition for the phase diagram of the attractive Fermi gas at finite imbalance. 

PACS numbers: 03.75.Ss, 03.75.Hh 
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I. INTRODUCTION 



The physics of single particles immersed in an environment 
is ubiquitous in physics. It appears, for example, in the large 
polaron problem where a single electron is dressed by its in- 
teraction with phonons [1] or in models for dissipation and 
decoherence in quantum mechanics [2, 3]. In recent years, 
new directions for exploring quantum many-body problems 
have been opened through ultracold atoms [4]. In particu- 
lar, for degenerate Fermi gases, the interaction strength can be 
tuned over a wide range using Feshbach resonances. This al- 
lows to study impurity problems in a fermionic environment. 
A specific example is a gas of fermionic ''Li, where the two 
lowest hyperfine states are populated in a highly imbalanced 
situation. For this system, recent experiments have shown that 
the minority atoms ('down spins') apparently form a liquid 
of quasiparticles [5]. Due to the strong attractive interaction 
to the up-spin Fermi sea, the associated quasiparticle weight, 
as determined from a sharp peak in the rf-spectrum, is found 
to vanish beyond a critical interaction strength. This transi- 
tion may be interpreted as one, in which a single | -Fermion 
immersed in sea of t-Fermions can no longer propagate as 
a quasiparticle but forms a many-body bound state with the 
Fermi sea. The existence of such a transition has been pre- 
dicted by Prokof'ev and Svistunov [6, 7]. Using a novel di- 
agrammatic Monte Carlo method, they have shown that for 
strong attractive interactions, a molecular state is energetically 
favored compared to one in which the single down-spin forms 
a polaronic quasiparticle in the up-spin Fermi sea. In the 
present work, we analyze this problem by a simple variational 
wavefunction. It provides an analytically tractable model for 
the physics on the molecular side, thus complementing the 
variational description put forward by Chevy [8] for the po- 
laronic quasiparticle. Our wavefunction gives a ground state 
energy that matches perfectly the results of the diagrammatic 
Monte Carlo method. Moreover, it describes correctly the 
three-body physics of repulsive atom-dimer interactions in the 
deep molecular limit and has zero residue for the down-spin 
Green function. The variational wave function is also used to 
determine the saturation field hg beyond which a two com- 



ponent Fermi gas is fully polarized and the behavior of the 
so-called contact coefficient introduced by Tan [9] in the limit 
of strong imbalance. 



II. THE FERMI POLARON AND ITS QUASIPARTICLE 
WEIGHT 

A simple variational wavefunction for the (NH-l)-particle 
problem of a single down-spin Fermion immersed in a sea 
of spin-up Fermions has been introduced by Chevy [8]. It is 
based on an expansion up to single particle-hole excitations 
around the unperturbed Fermi sea 



kq "q-k"k^q 



k,q 



\FS^) 
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Here and in the following sums on k and q with a prime are 
restricted to fc > kp and q < kp, respectively. Moreover, 
is the N-particle Fermi sea and the creation operators 
of up- and down-Fermions with momentum k are denoted by 
uj^ and (ij^. Despite the restriction to single particle-hole ex- 
citations, which is difficult to justify for the relevant case of 
zero range interactions that can create particle-hole pairs at 
arbitrary momentum, Monte Carlo calculations show that the 
ansatz (2. 1) gives a ground state energy that is very accurate, 
in particular at unitarity, where the scattering length a is infi- 
nite [6, 7]. The reason why the leading term in an expansion 
in the number of particle-hole excitations gives very good re- 
sults for the ground state energy can be traced back to the 
decoupling of higher order terms for vanishing hole momenta 
q = [10], i.e. contributions with more than one particle hole 
excitation interfere destructively. 

The wavefunction (2.1) describes the added down-spin as a 
quasiparticle dressed by its interaction with the up-spin Fermi 
sea. The virtual cloud of particle-hole excitations leads to a 
quasiparticle energy 



E{p) = Aep 
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(2.2) 



at low momenta |p| ^ kp that contains a 'binding energy' 



Aep < of a single down-spin to the Fermi sea and an ef- 
fective mass TO* [11]. Here, the Fermi energy is defined by 
Ef = kp/ {2m) (we use h ~ 1 throughout the paper) with 
a Fermi momentum kp that is related to the up-spin density 
by the standard relation = kp/{6Tr^) for a single compo- 
nent Fermi gas. Since we are interested in the limit of van- 
ishing down-spin density ni 0, these are the relevant en- 
ergy and momentum scales. The dimensionless coefficient A 
and the effective mass m* have been determined from varia- 
tional Monte Carlo calculations at the unitarity point [11] and 
from a T-matrix approximation at arbitrary values of the di- 
mensionless interaction strength i' = l/(fci?a) [12]. Very re- 
cently they have also been measured experimentally, giving 
A « -0.64(7) [5] and m*/m = 1.17(10) [13] at unitarity, in 
rather good agreement with the theoretical predictions. 

From a many-body point of view, the criterion that a sin- 
gle added down-spin is indeed a proper quasiparticle can be 
expressed by defining the quasiparticle residue Zi from the 
long-time limit 



= lim |Gx(p = 0,t)| ^0 

t — 



(2.3) 



of the down-spin Green-function at zero momentum. Within 
the variational wavefunction (2. 1) this residue is simply given 
by the probability = |(/)oP that an added down-spin 
at momentum p = is not mixed with plane waves at 
nonzero momenta q — k 7^ through particle-hole excita- 
tions. The fact that the coefficient l^oP of the Chevy wave- 
function coincides with the quasiparticle weight can be de- 
rived formally by noting that the ansatz (2. 1) is equivalent to a 
non-selfconsistent T-matrix approach for the down-spin Green 
function, which sums the particle-particle ladder for the ver- 
tex part r(k,a;) [12]. It is then straightforward to see that 
|(/)oP = |1 ^ 9ijS|i^=o|^^ coincides with the standard def- 
inition of the quasiparticle weight via the energy derivative 
of the down-spin self energy Y,{p,uj) at zero frequency and 
momentum. The numerical value of Zi at unitarity v = 
is Zi{v = 0) ~ 0.78 within the Chevy ansatz. This is much 
larger than the experimentally observed value Z[ = 0.39(9) 
which is likely to be a lower bound, however [5]. Smaller 
values Zi {v = 0) = 0.47 of the quasiparticle weight at unitar- 
ity are found from a 1/iV-expansion of the attractive fermion 
problem at strong imbalance, which is equivalent to a non- 
selfconsistent T-matrix approximation with the bare chemical 
potential [14]. 

In Fig. 1 we show the quasiparticle residue Zi for the 
minority Fermion as a function of = l/{kpa) within the 
ansatz (2.1) in comparison with the recent experimental re- 
sults [5]. Apparently, the expansion up to single particle- 
hole excitations considerably overestimates the quasiparticle 
residue even though it gives reliable results for the ground 
state energy. A much more basic shortcoming of the ansatz 
(2.1), however, appears if one considers the BEC-limit v ^ 1. 
Indeed, the ansatz predicts a finite value of Zi at arbitrary in- 
teraction strengths, even in the deep molecular limit. In this 
limit, however, an added down-spin will form a bound state 
with one of the up-spin fermions and can no longer propa- 
gate as a coherent quasiparticle. One thus expects that Z^ 
vanishes identically beyond a critical strength vm > of the 
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FIG. 1: (Color online) Quasiparticle residue of the minority 
Fermion as function of (kpa)^^, calculated using Chevy's varia- 
tional ansatz (2. 1). In the regime where the ansatz (2. 1) breaks down, 
Zi is drawn as dotted line. The red dots correspond to the experi- 
mentally measured quasiparticle residue from the MIT group [5] at a 
minority concentration of 5%. 



interaction, consistent with the experimental findings [5]. It 
is important to note, that the formation of the bound state is a 
genuine many-body effect at any finite density of the up-spin 
Fermi sea. Indeed, the binding is of a two-body nature only in 
the trivial limit v ^ I (that is, effectively, for kp ^ 0), where 
the bound state is formed with a single up-spin Fermion. By 
contrast, just beyond the critical value vm, the down-spin 
is effectively compensated by forming a singlet with many 
up-spin Fermions, somewhat similar to the physics of a lo- 
calized Kondo-spin interacting antifeiTomagnetically with a 
sea of conduction electrons at temperatures much below the 
Kondo-temperature [15]. Note, however, that in the Kondo 
problem the impurity spin is not fixed and the transition from 
an uncompensated spin to an effective singlet state appears as 
a continuous crossover from high to low temperatures. In the 
present problem, instead, there is a discontinuous transition in 
the ground state as a function of the attractive coupling v. 

An indication that the variational wavefunction (2.1) is not 
applicable for strong attractive interactions is provided by 
considering the Thouless criterion for a superfluid instability 
in which up-and down-spins are paired in an s-wave super- 
fluid [16]. Evaluating the relevant vertex function within the 
T-matrix approximation for arbitrary values of the down-spin 
chemical potential it is found that the Thouless criterion 
r^^(k = 0,cj = 0) = leads to a critical value for fi^ 
that is below the value /ij^ = ^^(p = 0) obtained from the 
ground state energy of the variational state (2. 1) provided that 
V > 1.27. A second argument that indicates the breakdown 
of the ansatz (2.1) in the regime w 3> 1 is the behavior of 
the ground state energy. Indeed, in a systematic expansion in 
powers of the scattering length a ^ 0+ the ground state en- 
ergy of a single added down-spin relative to the free up-spin 
Fermi sea is expected to be of the form 



E' 



(2.4) 
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Its leading contribution is just the molecular binding energy 
Ef, = —l/{ma^) < 0. The contribution ep of order a'^ 
accounts for the removal of one t-Fermion from the Fermi- 
sea that is required for the formation of the molecule. The 
last term, of order a is the mean field repulsion between the 
molecule and the Fermi-sea. Its interaction strength gad = 
Snaad/nT' is related to the exact atom-dimer scattering length 
aad = 1.18a that has first been calculated in connection 
with neutron-deuteron scattering [17] (for a recent derivation 
in the cold gas context see Petrov et. al. [18]). It turns 
out, however, that the variational ansatz Eq. (2.1) leads to 

E Eb - ep/^ + 0{a) which is too high by £^72 

compared to the exact asymptotics (2.4). The reason for this 
discrepancy can be seen easily from the structure of Chevy's 
wavefunction. On the BEC side, the dominant contribution 
comes from the q = terms, i.e. 

^'0k,odLkiiI^o|i^5f) (2.5) 

k 

describing the molecule formation of an up- and a down-spin 
with opposite momenta. This contribution is not optimal how- 
ever, since it creates a hole in the center of the f-Fermi-sphere. 
Energetically, it would be favorable to replace uo\FS^) with 
a (N-l)-particle Fermi sea |FS'j^~^), leading to a ground state 
energy that is lower by £i?. Within the ansatz (2.1), this would 
require terms with an arbitrary number of particle-hole exci- 
tations in order to reshuffle the Fermi-sea in such a way that 
the hole vanishes. 



III. VARIATIONAL ANSATZ IN THE MOLECULAR 
REGIME 

In order to describe the physics of bound state formation 
in the regime w ^ 1, we propose a variational ansatz for the 
(NH-l)-body problem that complements the ansatz (2.1) de- 
scribing a Fermi polaron with a finite quasiparticle residue. 
Our ansatz gives the exact behavior (2.4) of the ground state 
energy in the BEC-limit up to linear order in a. The associated 
variational wavefunction 

k 

k',k,q 

is a natural generaUzation of the Chevy ansatz and is con- 
structed by adding a (|, i)-pair to a {N — l)-particle Fermi 
sea of l-Fermions, together with the leading term in an ex- 
pansion in particle-hole excitations. Again, sums on k, k' and 
q are restricted to k,k' > kp and q < kp, respectively. 
The first term accounts for the formation of the molecule in 
the presence of the t-Fermi sea and gives the correct next-to- 
leading-order ground state energy in the BEC-limit, avoiding 
the problem of creating a hole in the t-Fermi sea. The single 
particle-hole excitation in the second term describes the lead- 
ing order contribution to the interaction of the dimer with the 



Fermi sea apart from Pauli-blocking effects, that are already 
accounted for in the first term. An important feature brought 
about by the inclusion of the second term in Eq. (3.1) is that 
it amounts to an exact treatment of the three-particle problem. 
Indeed, as is shown in detail in the Appendix, the set of cou- 
pled equations (3.6)-(3.9) that determine the coefficients of 
the variational many-body wavefunction reduce, in the three- 
particle limit, precisely to the integral equation for the exact 
solution of the three-body problem by Skorniakov and Ter- 
Martirosian [17]. As a result, the exact atom-dimer scattering 
length ttad = 1.18a appears in the asymptotic behavior of the 
ground state energy (2.4), giving rise to the correct next-to- 
next-to-leading order behavior of the ground state energy in 
the BEC-limit. 

Obviously, the ansatz (3.1) is not capable of describing the 
whole range of scattering lengths correctly. In particular, it 
does not capture the weak coupling limit a ^ 0^. Indeed, 
the |-Fermion in the first term is always added at momenta 
k > kp, leading to a ground state energy that is too high 
by ep in the weak coupling limit. Our ansatz (3.1) is there- 
fore complementary to the Chevy wavefunction (2.1), which 
correctly describes the situation at weak coupling up to and 
slightly beyond the unitarity limit. 

From a physical point of view, the two variational wave- 
functions (2.1) and (3.1) characterize very different ground 
states. Chevy's ansatz describes a Fermi polaron with a finite 
quasiparticle residue, which allows to build a normal Fermi 
liquid at a finite concentration of the down-spin Fermions, 
provided that interactions between the quasiparticles have no 
attractive channels (see section IV. below). By contrast, the 
wavefunction (3.1) describes a bosonic molecule interacting 
with a Fermi sea. At a finite concentration ^ 0, the result- 
ing ground state is expected to be a superfluid, coexisting with 
unpaired up-spin Fermions. The critical coupling vm, where 
the ground state energies of the two variational wavefunctions 
intersect is thus expected to separate a normal fluid from a su- 
perfluid ground state of the attractive Fermi gas in the limit of 
very strong imbalance. 

The variational ansatz (3.1) is based on a single channel 
model that describes the attractive interactions between both 
spin states. For computational purposes, however, it turns out 
to be easier to start from the more general two-channel model, 
which is defined by the Hamiltonian 

+-% E (^pCp_p' jcp,,i + /i.e.) . (3.2) 
p'.p 

Here, denotes the bosonic creation operator of a molecule 
with momentum p and cj, ^ are the fermionic creation oper- 
ators for the two species a —1^, j. The free particle disper- 
sion is denoted by ep = /{2m) and the factor 1/2 in the 
first term accounts for the factor two in the molecule to single 
Fermion mass ratio. The bare values of the detuning j/q and 
the Feshbach coupling strength go can be related to the physi- 
cal s-wave scattering length a and the interaction range r'o via 
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(3.3) 
(3.4) 



The two-channel Hamiltonian (3.2) is equivalent to a single 
channel model in the interesting limit of zero-range interac- 
tions ro ^ (i.e. for broad Feshbach resonances), as can be 
seen easily by integrating out the bosonic degrees of freedom. 

The corresponding variational ansatz to (3.1) in the two- 
channel model has two additional terms (~ 770, r/kq) where 
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the closed-channel state is occupied 
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Calculating the expectation value {ijjo\H — E\ipo), taking the 
derivatives with respect to the infinite set of variational param- 
eters rjo, ^k, ?7kq, ^k'kq and setting them equal to zero leads 
to the following set of coupled equations 
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(3.6) 
(3.7) 

(3.8) 
(3.9) 



Note that the ground state energy E is measured with respect 
to the N-particle Fermi sea, which explains the occurrence of 
the ep terms in the above equations. Moreover, using the N- 
particle Fermi sea as reference scale, the ground state energy 
E is equivalent to the chemical potential fi^ = E of the single 
down-spin. 



A. No particle-hole excitation 

Neglecting for a moment the contribution of particle-hole 
excitations in (3.5), i.e. setting rj^q = ^k'kq — 0, the ground- 
state energy is determined by the equations (3.6) and (3.7) 
alone. Performing the integrations and taking the zero-range 
limit ro 0, they reduce to a simple transcendental equation 



2kFa 
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(3.10) 



In the BEC-limit a 0+, Eq. (3.10) gives rise to a ground 
state energy of the form (2.4). The associated atom-dimer 
scattering length, however, is given by its value a^f " = 
(8/3) a in the Born approximation. More generally, it turns 
out that Eq. (3.10) is exactly equivalent to the Thouless- 
criterion F"^ (k = 0, cj = 0) = if the vertex function is cal- 
culated within a non-selfconsistent T-matrix approach where 
only the particle-particle ladder is summed, as discussed in 
section II. The resulting ground state energy is below that of 
the Fermi polaron if u > 1.27. 



B. Full variational treatment 

In the general case Tykq 7^ 0, ^k'kq 7^ 0, the equations 
(3.6)-(3.9) can be reduced to a single homogeneous Fredholm 
equation of the second kind for the variational parameters Tykq 
in the thermodynamic limit (again, the zero range limit has 
been taken already) 

Y,' K{E- k, q; k', q')r?k'q' = . (3.11) 

k%q' 

The associated Kernel K{E\ k, q; k', q') is given by 

y<5k,k' 1 ^^<5q,q' 
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(3.12) 

(3.13) 
Eq (3.14) 

(3.15) 
(3.16) 



Due to the isotropy of the system, the variational parameters 
77kq = ?/(^, 9, cos ^kq) depend only on the magnitudes of the 



two momenta k and q and the angle between them. This al- 
lows Eq. (3. 1 1) to be reduced to a three dimensional integral 
equation. 

The ground state energy E is now simply obtained by the 
condition that the Fredholm determinant of the kernel K van- 
ishes. We evaluate the Fredholm determinant numerically 
by discretizing the integral equation using a Gauss-Legendre 
quadrature and calculating the determinant of the correspond- 
ing linear equation system. The order of the quadrature for the 
k, q and cos6'kq integral were chosen as 11, 11 and 4, lead- 
ing to an error of 10^'' of the ground state energy at the 
unitarity point a oo, where the convergence is slowest. 

The ground state energy as function of (kpa)^^ is shown 
in Fig. 2. Apparently, our ansatz (3.1) leads to a ground 
state energy that is below that of Chevy's ansatz for interac- 
tion strengths larger than (kpo)'^^ = 0.84. This value is in 
good agreement with the diagrammatic Monte-Carlo results 
by Prokof'ev and Svistunov [6, 7], who obtain (kpa)^^ = 
0.90. In fact, the small discrepancy is entirely due to the fact 
that we intersect our molecular ground state energy with that 
obtained using the Chevy wavefunction, which is not precise 
near vm- The Monte Carlo results in turn give better values 
for the polaron energy, shifting the intersection slightly to- 
wards the BEC regime, as can be seen in Fig. 2.. Yet, as far 
as the molecular ground state energy is concerned, our results 
agree perfectly with the Monte Carlo data, down to the small- 
est coupling V ~ 0.6 where they have been calculated. 

It is interesting to note, that the approximation q = in the 
wavefunction (3.5) (i.e. pinning the hole-wavevector at zero 
momentum), leads to a ground state energy that differs from 
the calculation with the full wavefunction by at most 3% in the 
regime (fci?a)^^ > 0.84, where the ansatz is valid. The sit- 
uation thus appears similar to that in the polaron case, where 
Combescot et al. [10] have shown that an expansion in hole- 
wavevectors works very well for Chevy's ansatz at unitarity. 



C. Quasiparticle Residue 

We now show that the quasiparticle residue Z| of the 
Fermion, which can be thought of as a kind of order parameter 
of the transition from the polaron to the molecular state, van- 
ishes identically in the thermodynamic limit for the variational 
wavefunction (3.5), that gives a lower ground state energy on 
the molecular side of the critical coupling v m ■ 

Since the variational ground state wavefunction does not 
allow to calculate the full down-spin Green Function, the def- 
inition (2.3) of the quasiparticle residue is not applicable. In- 
stead, we use the standard connection between Zi and the 
jump in the momentum distribution at the Fermi momentum 
kpi, and the latter is zero in the limit of a single down-spin. 
The momentum distribution of the j-Fermion within the vari- 
ational ansatz (3.5) is given by 

"i> = I^pI" + 2 l^k'kql' 5p,q-k'-k (3.17) 

k'kq 
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FIG. 2: (Color online) Ground state energy E — Et (binding en- 
ergy El, — —l/{ma?) subtracted) in units of the Fermi energy ef 
as function of {kFa)~^. Blue solid line: Chevy's ansatz (2.1); red 
line with full squares: ansatz Eq. (3.1); black dashed line: BEC- 
asymptotics (2.4); orange dash-dotted line: Thouless criterion (3.10). 
The open black diamonds and green triangles correspond to the QMC 
results for the molecule- and the polaron energy from Prokof'ev and 
Svistunov [6]. 

and is normaUzed via 

i = E"p = E'i^kiV2^'ia'kqi' . (3.18) 

p k k'kq 

The normaUzation condition requires the coefficients to scale 
with the system volume as ~ 1/ \/V and ^k'kq ^ 1 /V^^"^. 
Since an upper bound to the quasiparticle residue Z| is given 
by the momentum distribution at p = and = for p < 
kp,we find that 

< nl^, = 2^' l^'kql' <5q,k'+k ~ . (3.19) 

k'kq 

As a result, the quasiparticle residue of the molecular 
wavefunction scales inversely with the volume of the system 
and thus vanishes in the thermodynamic limit. This is in con- 
trast to Chevy's wavefunction, where Z^ = |(/)oP is always 
finite. The two wavefunctions (2.1) and (3.1) therefore indeed 
describe qualitatively different ground states. In particular, no 
sharp peak is expected in the minority rf-spectrum at coupling 
strengths v > vm, consistent with the experimental observa- 
tion [5]. 

In the q = approximation, which captures the essen- 
tial properties of the variational ansatz (3.1), the quasiparticle 
residue Z| in fact vanishes identically. Indeed, 

<2^'|e-kko|' = 0, (3.20) 

k 

since, as can be seen from Eq. (3.9), the coefficients ^_kko 
f^ko — ?/-ko = vanish because r/ko only depends on the length 
ofk. 



IV. CONTACT COEFFICIENT AND PHASE DIAGRAM 



The analysis of the polaron to molecule transition in the 
previous section leaves two important questions open: what is 
the nature of the transition and what are its implications for 
the phase diagram of the strongly imbalanced gas? Now for 
the case of a single down-spin in an up-spin Fermi sea, the 
transition from a polaronic to a molecular state is a first or- 
der transition, where the quasiparticle residue exhibits a 
discontinuous jump from a finite value to zero at the critical 
coupling vm ~ 0.9. This is a result of the fact that the ener- 
gies of the two ground states, which have different quantum 
numbers, cross with a finite slope at vu (see Fig. 2.). It is im- 
portant to note that this crossing is not an artefact of extending 
the different variational states beyond their domain of validity. 
Indeed, as shown by Prokof'ev and Svistunov [6, 7], both the 
polaronic and the molecular state exist as stable excitations 
for V > Vm or v < respectively because the phase space 
for decay vanishes linearly with the magnitude of the energy 
difference. Both states are thus reachable as metastable con- 
figurations coming from the weak coupling or the molecular 
side, as expected for a first order transition. 

A different perspective on the first order nature of the po- 
laron to molecule transition is provided by considering the so 
called contact coefficient C. As shown by Tan [9], the mo- 
mentum distribution of Fermi gases with zero range interac- 
tions generically decays with a power law na{k) C/k'^ for 
large momenta. The associated coefficient C is identical for 
both spin components a =ti i [20] and is a measure of the 
probability that two Fermions with opposite spin are close to 
each other [21]. Using the adiabatic theorem derived by Tan 
[22] , the contact density can be determined from the derivative 
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of the ground state energy density u = E/V with respect to 
the inverse scattering length. Now the definition of the down- 
spin chemical potential yUj implies that the energy density u 
of the strongly imbalanced Fermi gas ^ 7i| to linear order 
in the minority density 7i| is of the form 



u = -Eft"-! + ^J■l''H + 



(4.2) 



where the first term is simply the energy of a non-interacting 
gas of spin-up Fermions. The dimensionless contact coeffi- 
cient s defined hy C = s ■ kpkp^ for a strongly imbalanced 
Fermi gas can thus be obtained from the derivative 



s = 



1 di-ni/ep) 



dv 



(4.3) 



of the negative down-spin chemical potential in units of the 
Fermi energy with respect to the coupling constant v. Since 
fil is precisely the energy E associated with adding a sin- 
gle down-spin, our result for the ground state energy of the 
{N + l)-particle problem immediately gives the contact den- 
sity of an almost fully polarized attractive Fermi gas (note that 




0.25- 



1.5 



FIG. 3: (Color online) Dimensionless contact coefficient s as func- 
tion of (/cFfl)^^, calculated from the two variational wavefunctions 
(2.1) and (3.1). Within this approach s is discontinuous at the critical 
coupling Vm. The black dashed line marks the asymptotics in the 
molecular limit, where s — 4v/3n is fixed by the two-particle bound 
state wavefunction in momentum space. 



this applies even on the molecular side v > vm, where the sin- 
gle added down-spin is not a propagating quasiparticle). The 
associated dimensionless constant s is shown in Fig. 3. It in- 
creases monotonically from weak coupling to unitarity and up 
to the critical coupling vm- At this point, there is a discon- 
tinuous jump upwards, that reflects the transition to a molec- 



ular state. Note, that the proportionality C ^ kp^ 



of 



the contact to the down-spin density makes C vanish in the 
limit of full polarization. This is expected, because the fully 
polarized gas at zero temperature is an ideal Fermi gas, with 
no tails in the momentum distribution. Apart from the jump 
at Vm, the behavior of the dimensionless contact coefficient 
s is rather c^ose to that obtained for the contact coefficient 
C = s ■ kp of the balanced superfluid along the BCS-BEC 
crossover [23] (note that the Fermi momentum fc^ of the bal- 
anced gas is related to that of the up-spin component used here 

~ 3 

by kp = hp (1 + <t), where a ~ (rif — ni)l{n^ + nj^) is the 
degree of polarization at a fixed total number of particles). In- 
deed, in weak coupling one obtains .s„c = (2/37rv)^ from the 
mean field attraction of the polaron to the up-spin Fermi sea, 
while s ~ 0.08 at unitarity and sbeg = '^v/'in in the molec- 
ular limit, very similar to the behavior that is found for s in 
the balanced superfluid [23]. 

The solution of the (A^ + l)-body problem for arbitrary cou- 
pling strengths v has also implications for the phase diagram 
of the imbalanced Fermi gas in the regime near complete po- 
larization. For a discussion of this issue, it is convenient to 
introduce an effective magnetic field h that couples to the two- 
different spin-states a i in the standard form 



H' 



h I I - iV, 



(4.4) 



of a 'Zeeman' field that favors a finite population imbalance 
a ~ (n| — 7i^)/(n-f + nj^) > 0. At a fixed total density n, the 
ground state energy u per volume is then a function of n and 
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M. 



SF, 



molecular side, Eq. (2.4) gives 



-0.5 



0.5 



(kpa) 



FIG. 4: Qualitative pliase diagram of tiie imbalanced Fermi gas as a 
function of tiie inverse coupling strength (kpa)^^ and the effective 
magnetic field h/ep- The thick line indicates a first order phase tran- 
sition and the different phases are labeled as in [25], i.e. Nf^: fully 
polarized normal phase, A'^ppi partially polarized normal phase, SFq: 
balanced superfluid, SF-p : polarized superfluid. The points M and 5* 
are discussed in the text. The precise structure of the phase diagram 
in the nontrivial regime he < h < hs is likely to contain unconven- 
tional superfluid phases in addition to the TVpp and SFp phase, which 
are not shown in our figure. 



h. It determines the chemical potentials of the majority and 
minority species from = fi± h where yu = du{n, h)/dn 
is the average chemical potential. In addition, it also fixes the 
imbalance from — = —du{n,h)/dh. The choice of 
an ensemble with fixed values of n and h is convenient for a 
discussion of the ground state phase diagram of the attractive 
Fermi gas at arbitrary coupling v, both in the homogeneous 
case and in the presence of a harmonic trap [24]. Indeed, there 
are two critical fields hc{v) and hs{v) that separate two sim- 
ple limiting phases from a regime, in which nontrivial ground 
states are expected: The lower critical field he is defined by 
a{h) = for h < he and determines the boundary of the 
balanced superfluid phase (denoted by SF^ in Fig. 4, follow- 
ing the notation used by Pilati and Giorgini [25]). The upper 
critical (or 'saturation') field hs, in turn, is defined by the con- 
dition of complete polarization cr(/i) = 1 for h > hg. Since 
a single component Fermi system has vanishing interactions 
in the ultracold limit, this regime is just an ideal Fermi gas, 

i. e. it is a normal, fully polarized state. The qualitative struc- 
ture of the zero temperature phase diagram as a function of 
the interaction parameter v = 1/ {kpa) and the effective mag- 
netic field h in units of the bare Fermi energy e^^ of the fully 
polarized gas is shown in Fig. 4. 

In this diagram, the upper line hs{v) is completely fixed by 
our calculation above of the energy /i| associated with adding 
a single down-spin to an up-spin Fermi sea. Indeed, since 
/i| = £i? along this line, we have hs = (ej- — ^|)/2. In 
terms of the constant A{v) introduced in Eq. (2.2), this leads 
to hs/ep = (1 — ^)/2, giving hs ~ O.Slei? at unitarity from 
the precise numerical value of the polaron energy [7]. On the 



1 



2'nv 



(4.5) 



which is very accurate even at u = vm- The point M along 
this line separates a regime where a single down-spin is a well 
defined fermionic quasiparticle from one, in which it is bound 
to the up-spin Fermi sea. The first order nature of the transi- 
tion shows up as a discontinuity of the slope in hs{v) at M 
which is, however, hardly visible in Fig. 4. For a finite den- 
sity of down-spins, the point M appears as an endpoint of a 
line that separates a phase with a finite Fermi surface volume 
1.5 ri| ^ to its left from one with fij^ = [26]. Using the 
generalized Luttinger theorem derived by Sachdev and Yang 
[27], the expected polarized superfluid (SFp) phase on the 
molecular side has a condensate of 'dimers' plus an up-spin 
Fermi sea, whose volume il| = (27r)'^(7i| — rt^) is set by the 
imbalance. This is consistent with the naive picture that the 
density of unpaired up-spins is simply n| — even though 
the 'dimers' in the vicinity of the transition are far from lo- 
cal (t, j)-pairs. In principle, this simple picture of the SFp- 
phase as a BEC coexisting with a sharp, single Fermi sur- 
face of unpaired up-spins is unstable with respect to p-wave 
pairing due to the induced interactions between the unpaired 
fermions through the superfluid [28]. In practice, the nontriv- 
ial superfluid phase of the unpaired up-spins is exponentially 
suppressed for strong imbalance. Moreover, quantitative re- 
sults for the p-wave instability can be derived only in second 
order in 1/u ^ 1, where the resulting energy scales are ex- 
ponentially small compared with ep- In practice, therefore, 
the phase with p-wave pairing among the unpaired up-spins 
seems hardly accessible experimentally. 

A nontrivial issue that has been neglected in the discus- 
sion so far is the question whether a gas of polarons or bound 
molecules is indeed stable at low but finite densities n i . On 



i- 

the weak coupling side, there is again an induced attractive in- 
teraction in the p-wave channel among both the up-spins and 
the down-spins, mediated by the other species. The ground 
state is thus expected to be a two-component p-wave super- 
fluid and not a normal Fermi liquid state [28]. Similar to 
the situation in the BEC-limit, however the energy scale for 
this instability is exponentially small in the regime where 
the calculation can be controlled. More importantly, as has 
been shown recently by Nishida [29], the effective interac- 
tion between two heavy down-spin fermions immersed in an 
up-spin Fermi sea is attractive in the p-wave channel only 
for weak coupling. Approaching unitarity, the p-wave inter- 
action becomes repulsive. Assuming that this result carries 
over to the relevant case of equal masses of the up- and down- 
spin Fermions, a finite density gas of down-spins will indeed 
form a normal Fermi liquid at unitarity, as was implicitly as- 
sumed in the calculations of the equation of state and den- 
sity profiles of the unitary gas beyond the critical imbalance 
(Tc — 0.4, where the balanced superfluid is no longer stable 
[11, 30]. On the molecular side, the phase immediately below 
the saturation field line hs{v) is expected to be a superfluid 
of (t, 4)-pairs at a very low density immersed in an 

up-spin Fermi sea. The fact that the atom-dimer repulsion 
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aad = 1.18 a is much larger than the dimer-dimer repulsion 
o-dd = 0.6 a [18], however indicates that a low density gas of 
molecules tends to phase separate from the up-spin Fermi gas. 
This phase separation has indeed been found from an extended 
BCS-description of the BCS-BEC crossover in an imbalanced 
gas [31, 32, 33, 34]. It has recently been seen also in the vari- 
ational Monte Carlo calculations by Pilati and Giorgini [25]. 
Their results indicate that a section between vn — 0.73 and 
a triple point at vt — 1.7 along the /ig-line is actually a first 
order line, where the polarized superfluid disappears with a fi- 
nite jump in density as the effective field h increases through 
hs- As shown above, the point M lies in the interval between 
Vn and vt and thus the polaron to molecule transition would 
not be accessible at any finite minority density, at least not in 
an equilibrium situation. Clearly, our variational calculation 
for the single down-spin problem cannot address the question 
of phase separation. An unexpected feature of the hg-line in 
the presence of phase separation is the fact that the transition 
across hg is predicted to be continuous up to v^, first order be- 
tween I'N and Vt and continuous again for v > vt- The rather 
large value v-r — 1.7 up to which phase separation is pre- 
dicted also appears surprising. Indeed, in this regime a mean 
field theory describing a Fermi gas coexisting with a EEC of 
molecules gives for the energy per volume as a functional of 
the density difference Sn = ri^ — ni and the dimensionless 
field h — h/ep the simple form 



u(ri, Sn, h) 



£f 



3 /(5n\2/3 
51-) ^" 



— Sn) 



CLad {n - Sn) Sn ddd {n - Sn)'^ 



2nv 

h Sn 



12nv 



(4.6) 



Here, Gad = ciad/a and add = o-dd/a are the atom-dimer and 
dimer-dimer scattering length measured in units of the atom- 
atom scattering length and v — {kpa)~^. The true ground 
state energy density u{n, h) as function of the total density n 
and the effective magnetic field h is determined by the min- 
imum of the Landau energy (4.6) with respect to the density 
imbalance 



diulep) /Sn\'^/^ o 
(-) + ^' 



aad _ ^Sn 



\ n / 



h = 



(4.7) 



d Sn \ n J 2'Kv \ n 

_l_ add fSn 
Qtiv V n 

This equation determines the imbalance cr = as a func- 
tion of the field h and indeed, it correctly describes the ex- 
act asymptotic results for both the saturation field hg and the 
lower critical field he in the limit v ':$> 1 (see Eqs. (4.5) and 
(4.9)). In this simple model, phase separation between a polar- 
ized superfluid phase and a fully polarized Fermi gas appears 
for coupling constants v < Vc,ps, below which the energy 
density (4.7) has a second minimum at full polarization a = 1. 
This occurs at 

^c.PS = ^(«a.-^). (4.8) 

With the exact values dad = 1.18 and ddd = 0.6 one ob- 
tains Vc.ps = 0.516, where the simple expansion (4.6), how- 



ever, is no longer valid. From the calculation above, the triple 
point vp beyond which phase separation appears in an almost 
fully polarized gas lies at a much smaller value of the cou- 
pling strengths than found previously [25, 34]. At finite tem- 
peratures, phase separation is suppressed by the presence of a 
mixing entropy, which may explain that it is not observed in 
the experiments, where T w 0.15Ti?. 

Concerning the lower critical field hc{v), its weak-coupling 
limit is determined by the well known Chandrashekar- 
Clogston result he ~ /S./ a/2, beyond which the balanced 
BCS-pairstate is unstable [35, 36]. In the molecular Umit 
V 1, the critical field 



he 
£f 



1 



2'Kva 



aad 



add \ 
6 / 



(4.9) 



follows essentially the two-particle binding energy with cor- 
rections due to the atom-dimer and dimer-dimer scattering 
lengths aad and add, respectively. At unitarity, he = 0.26ep 
is a universal constant times the bare up-spin Fermi energy 
ep [11]. As a result, there is a wide range hg/hc = 3.12 
between the balanced superfluid and the fully polarized gas, 
much larger than that found in an = oo theory of the im- 
balanced attractive Fermi gas, where hg/hc = 1.24 [37]. For 
h > he the balanced superfluid is destroyed by the onset of 
a finite polarization a ^ 0, which leads to a mismatch of 
the Fermi energies. An effective field theory due to Son and 
Stephanov [38] indicates, that the phase beyond the balanced 
superfluid exhibits a spatially varying superfluid order of the 
FFLO type as also found at weak coupling. The transition is 
first order with a jump both in polarization a an total density, 
a situation that is also found in the case of a direct transition 
between a balanced superfluid and a partially polarized nor- 
mal phase [11, 30]. In an ensemble with a given density that 
is used here, the /ic-line would then split into two distinct lines 
hci and hc2, as noted by Sheehy and Radzihovsky [32]. In our 
diagram in Fig. 4, he denotes the boundary of the balanced su- 
perfluid at given density, which is well defined without spec- 
ifying the state that is reached at nonzero polarization. The 
first order nature of the transition is found only up to a split- 
ting point S, beyond which the fermionic excitations have 
their minimum at p = 0. When this is the case, additional 
up-spin Fermions can be added by filling up a Fermi surface 
whose volume fl^ ^ a ^ (h — hc)^/^ increases continuously 
from zero. The transition from the balanced superfluid to a 
polarized superfluid with unpaired excess Fermions is there- 
fore continuous and preserves superfluidity. The precise lo- 
cation of the splitting point S has been determined recently 
from a calculation of the fermionic excitation spectrum along 
the BCS-BEC crossover of the balanced gas [23]. It is located 
at vs — 0.63 and hc{vs) = A ~ 0.6£f, at considerably 
larger coupling strengths than predicted by mean-field theory 
where the splitting point coincides with the zero crossing of 
the chemical potential (note the factor 2^/'^-difference with 
the result in Ref. [23], which is due to the fact that the up-spin 
Fermi wave vector and not that of the balanced case appears in 
our present coupling constant v). This is in agreement with the 
calculation of the splitting point within an e — 4 — d-expansion 
by Nishida and Son [39] but is larger than the value vs — 0.5 
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found for the splitting point in the Monte Carlo calculations of 
Pilati and Giorgini [25]. Note that the possibility of extracting 
the critical coupling of the splitting point from a calcu- 
lation of the balanced gas relies on the fact that the ground 
state energy is independent of the field h in the whole regime 
h < he because the polarization a = -~ndu{n, h)/dh van- 
ishes. The nature of the phase diagram near the splitting point 
has been discussed by Son and Stephanov [38] using an effec- 
tive field theory. In particular, the phase immediately beyond 
he is expected to be of the FFLO-type, with a spatially os- 
cillating superfluid order parameter that appears also in weak 
coupling beyond the Chandrasekhar-Clogston limit [40]. It is 
an open question, how this nontrivial superfluid evolves into 
a normal phase in which the two spin components each form 
a Fermi liquid. In fact it is this latter phase, which describes 
the experimentally observed density profiles [41] at unitary 
extremely well [11, 30]. It is also an open issue, of how to 
separate in detail the regime between the lower critical field 
and the saturation field into a regime where an imbalanced 
Fermi liquid or a polarized superfluid phase appear as ground 
states. In the phase diagram of Pliati and Giorgini, the first or- 
der line that bounds the balanced superfluid up to the splitting 
point S extends as a first order line up to hg at the coupling 
vn and then continues along hg up to the tricritical point vt- 



V. CONCLUSIONS 

From a variational wavefunction that describes the (A^+ 1)- 
particle problem of a single down-spin interacting strongly 
with an up-spin Fermi sea, we have discussed the physics of 
the strongly imbalanced Fermi gas. In particular, we have fo- 
cused our attention on the quasiparticle residue and the con- 
tact coefficient C. The latter exhibits a discontinuous jump 
at the polaron to molecule transition which might be detected 
by measurements of the closed channel fraction similar to the 
analysis of the experiments by Partridge et.al. [42] due to 
Werner, Tarruell and Castin [43]. A motivation for this work 
were the recent experiments of Schirotzek et.al. [5], who have 
observed a transition from a Fermi liquid phase of polaronic 
quasiparticles near unitarity to a phase in which the quasi- 
particle residue vanishes. As shown above, this transition is 
expected to be a discontinuous one in the single down-spin 
limit. Apparently, however, vanishes in a continuous man- 
ner in the experiment (see Fig. 1 ). Apart from the uncertainty 
in extracting Z| from the sharp structure in the minority rf- 
spectrum, this discrepancy is probably due to the fact that the 
Chevy wavefunction (2.1) strongly overestimates the quasi- 
particle residue near the polaron to molecule transition. A 
reliable quantitative calculation of the rf-spectra for both finite 
concentrations and at finite temperatures is, unfortunately, not 
available. The existence of a stable finite density gas of po- 
larons in the regime up to w ~ 0.9 however indicates that the 
interaction between them is repulsive, so that they indeed form 
a Landau Fermi liquid below the critical coupling vm- As dis- 
cussed in section IV, the detailed structure of how this phase 



connects to the nontrivial superfluid phases expected near the 
splitting point S and on the BEC-side is a major and still open 
problem. 
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APPENDIX: THREE PARTICLE LIMIT 

In the following we briefly show, how the exact solution of 
the three Fermion problem can be obtained from the integral 
equation (3.11) for the full variational wavefunction (3.5). If 
the (N-l)-particle Fermi sea is reduced to a single |-Fermion, 
only the q = terms remain in the variational wavefunction 
(3.5). Thus, starting from the integral equation (3.11) in the 
thermodynamic limit, taking the limit fc^ and setting 
q = 0, one arrives at the simplified equation 



ako '7ko 



Vk'O 



(2^)3 Sk'kO 



(A.l) 



Inserting the coefficients ako and i?k'ko from Eqs. (3.14) and 
(3.15) explicitly, the integral equation (A.l) takes the form 



3fc2 



- mE 



(27r)3 k^ + k'^ + - mE 



(A.2) 



Note that 7]k = ?/ko corresponds to the Fourier transform 
of the relative wavefunction between the (|, J,)-molecule and 
the additional f-Fermion. The integral equation (A.2) is ex- 
actly the same as the one obtained by Skorniakov and Ter- 
Martirosian [17] for the three Fermion problem. In particu- 
lar it is equivalent to Equ. (29) in [17], which corresponds to 
three-nucleon scattering with total isospin T = 1 /2 and to- 
tal spin S — 3/2 (note that the spin part only contributes an 
unimportant prefactor to the wavefunction in this case). 
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